A mechanism for phase separation in copper oxide superconductors 
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A two-band Hubbard model is used to describe the band structure and phase separation (PS) 
in multiband superconductors, especially in cuprates. We predict a large peak in the density of 
states at the Fermi level in the case of optimum doping, corresponding to the minimum energy 
difference between the centers of two hole bands. For strong interband hybridization, a metal- 
insulator transition occurs near this optimum doping level. We suggest a mechanism of PS related 
to the redistribution of holes between two Hubbard bands rather than to the usual antiferromagnetic 
correlations. We show that the critical superconducting temperature Tc can be about its maximum 
value within a wide range of doping levels due to PS. 

PACS numbers: 74.72.-h,71.27.+a,64.75.+g 



Nanoscale spatial variations in the electronic charac- 
teristics of different high-Tc superconductors have been 
commonly observed in the form of stripes gran- 
ular droplet-like structures four- unit-cell periodic 
(checkerboard) patterns or even more intricate ar- 
rangements y. Several theoretical models also sug- 
gest inhomogeneous electron structure in superconduct- 
ing cuprates and related materials 0, @, Moreover, 
the phenomenon of self-organized electron inhomogene- 
ity or phase separation (PS) is common to many strongly 
correlated electron systems jli,i9|]. For example, PS in the 
form of droplets, stripes, and checkerboard patterns has 
been observed in doped magnetic oxides (manganites) [ij . 
Many theoretical approaches [J |3j B @| explain PS as a 
result of the competition between electron localization, 
due to antiferromagnetic correlations, and delocalization 
in the non-magnetic (or ferromagnetic) regions. 



Here we consider an alternative mechanism (without 
antiferromagnetic correlations) for the PS arising in a 
Hubbard model with two (or more) bands. Even though 
this model is commonly used to describe high- Tc super- 
conductors PS has not been studied using such an 
approach. The two-band Hubbard model allows us to de- 
scribe PS of the droplet type (including the droplet size) 
observed in experiments (see, e.g., 2, 4]). In contrast to 
the usual mechanism, based on antiferromagnetic corre- 
lations, here we show that our proposed mechanism of 
PS gives rise to a density of states (DOS) at the Fermi 
level, which corresponds to the optimum doping in one 
of the phases. This results in a weak variance of the 
critical temperature within a broad doping range, as ob- 
served in experiments. The number of charge carriers in 
superconducting cuprates is not directly determine d by 
the doping level and can depend on temperature l3, U | ■ 
This fundamental problem can be also understood in the 
framework of two-band models Uj . 



Two-band Hubbard model. — The usual model used to 
describe superconducting cuprates is based on the Hub- 
bard Hamiltonian with three bands: Cu(3d), 0(2^2;), and 
0{2py) [ic|. In earlier studies, this Hamiltonian was re- 
duced to an effective single-band Hubbard Hamiltonian 
(see, e. g . . fiol fl^ ). However, recent computations (see, 



e.g. 
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) show that two distinguishable Cu-0 bond 
lengths in the CuOg octahedron, and the direct tun- 
nelling of holes between the oxygen atoms lead to the 
effective two-band Hubbard Hamiltonian jlij : 
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Here, ajj^.^. and Onaa are the creation and annihilation 
operators for holes in the state a — {p, d} at site n 
with spin projection a (a, a denote not-a and not-cr), 
the symbol (...) denotes a summation over the nearest 
sites, fi is the chemical potential, e^j is the energy dif- 
ference between the centers of the d and p bands and 
Cp = 0. The first term in Eq. ([!]) is the kinetic energy 
of the conduction holes, the second term is due to the 
chemical potential and the shift between the centers of 
d and p bands. The last two terms correspond, respec- 
tively, to the intra- and interband on-site Coulomb re- 
pulsion. In agreement to Ref. [ij, we assume that the 
Coulomb interaction is strong enough, that is C/", U' ^ 
t"^ , ed- Applying the mean field (so-called Hub- 
bard I) approximation, {famaait)nnpa'{t)a'laacr(to)) 
(nn/3CT')(2^«ma(T(0«noaa(*o)) (T is timc-orderiug opera- 
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tor), for Hamiltonian we derive the relationship 



e") G„;3(a;, k) =g,|l + ^£"'^(k)G^^(a>, k)J 
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for the one-particle Green's functions Gafj,aa' (n — 
no, t - to) = -~i{fanaa{t)al^ij^,{to)) in the frequency- 
momentum (cj, k) representation; where ffa = 1 — 
ria — ".a/2. The form of the function £:"^(k) de- 
pends on the symmetry of the crystal lattice. We an- 
alyze here a simple cubic lattice. In this case, we ob- 
tain e"'^(k) = ^"'^((k), w"'' = zi"^, and C(k) = 
— [cos{kid) + cos(fc2<i) + cos(fc3(i)] /3, where d is the lat- 
tice constant. Below we consider a purely paramagnetic 
state, that is, Ua] = Uai = and neglect the values 

(ajjag). Note, however, that the latter assumption does 
not affect much the obtained results [15| . 

Relationship ^ is the linear set of the equations for 
Gap, which can be easily solved. The calculated Green's 
functions determine the DOS and the energy of the 
system. Following this approach, we derive the DOS, 

in the form 

d^k 
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where ej{({'kj) and z;^(C(k)) are the eigenvalues and 
eigenvectors of the matrix e"'^ (^) — ^gagp w^^C—Sape" ■ 
Solving this eigenvalue problem, we obtain the energy 
spectrum Sj{li) of holes in two new bands (labelled by 
j = ±1) and coefficients determining the transforma- 
tion from p and d holes to the quasiparticles in these 
bands 
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jJ[{w<>'°- - iD'''')C + e]2 +4(w«'')2 \ , (4) 
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where Hj"^ — y/gagpw"^ . In contrast to p and d holes 
with short lifetime due to interband transitions, new 
quasiparticles with spectrum (j4|) have a longer lifetime 
and are scattered by, e.g., phonons and impurities. The 
transformation performed above is similar to the well- 
known Zhang-Rice derivation of the single-band effec- 
tive Hamiltonian if tpp = We denote the lower 
band as j = 1 and the upper one as j = —1. We 
can write paf3 using the dimensionless DOS po{E') = 
J S{E' — (^(k))c?'^k/(27r)'^ of uncorrelated electrons, 



Po(C)l 
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where Q (E) is the inverse function of Sj (E) . The number 
of electrons in the state a is 



2 / dE paa{E,na,n{) , 



(7) 



where ^-axin = £i(C = — 1) (note that pap depends on tIq, 
through the functions go)- The total number of charge 
carriers per site is n = X^q"-"!/^)- This equality, along 
with Eqs. ([1])-([7|), form a set of equations for the calcu- 
lation of Ua and /i. Performing a transformation of the 
operators Oq, employing the coefficients ([S]), we can find 
the Green's functions corresponding to the band spec- 
trum (lU and derive the DOS, pj for two j-bands 
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Band structure and density of states. — Superconduct- 
ing cuprates vary from strongly anisotropic to nearly 
isotropic materials. To account for this, we consider 
two limiting cases: two-dimensional square and 3D cubic 
symmetries. The obtained results are quite similar and 
below we present the band structure and DOS (Figs. 1 
and 2) calculated using Eqs. (|4]) and ^ only for the case 
of simple cubic lattice. Two very different quasiparti- 
cle spectra are shown in the insets of Figs. la,b. The 
anticrossing of two bands shown in the inset in Fig. la 
corresponds to a metallic behavior for any doping. When 
doping increases, the chemical potential p (shown by the 
dotted green line) shifts upward: at low doping we have 
one metallic band and one empty band; then two metal- 
lic bands; by further increasing the doping, one metallic 
and one filled band. For larger interband hybridization, 
we obtain a transition to insulator at some 
doping level. Indeed, for some doping, p, is located in the 
gap between zones. Nearby the anticrossing point of the 
two bands, the spectrum £j(C(k)) becomes flatter. This 
could indicate the so-called nesting of the Fermi surface. 
For energies close to anticrossing points, the DOS ex- 
hibit peaks (Fig. 1), which are large, pj cx I/^/Eq — E, 
when n is close to 1 in the vicinity of the metal-insulator 
transition (Fig. lb). The optimal doping for supercon- 
ductivity corresponds to the case when p is close to the 
energy where the peak of the DOS is observed. 

Phase separation in cuprates. — The energy, 
Shorn (n) = 2J2j'i^JEEpjiE), of the homoge- 
neous state, versus doping n, is shown in the inset of 
Fig. 2a. The curvature of £^hom('^) is negative between 
the two marked points ni and n2. This indicates the 
instability of the homogeneous state with respect to 
the separation into two phases with hole densities ni 
and 712, ni < n < n2 (see, e.g., Ref. T^. The energy 
Eps{n) = cEhom{ni) -I- (1 - c)Ehomin2) of the PS 
state with relative phase fractions c and 1 — c, which 
are determined by the charge conservation condition 
n = cni -|- (1 — c)n2, is lower than £'hom("-) between ni 
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FIG. 1: (Color online) Density of states p±\ versus energy 
E for quasiparticles (holes) located in bands j = ±1 for two 
different types of spectra, shown in the corresponding insets. 
The gapless spectrum shown in the inset (a) was calculated 
at n = 0.6, Ed ~ 0.15 eV, Wpp = 1 eV, Wdd = 0.1 eV, 
and relatively small interband hybridization Wpd ~ 0.25 eV. 
The gapped spectrum in the inset (b) was calculated for 
the same parameters but for stronger interband hybridiza- 
tion Wpd = 1 eV. In case (b), the transition to the insulating 
state occurs at a certain doping level for which the chemical 
potential /i (shown by the green dotted line) is located inside 
the gap. In both large peak in the DOS is observed 

at energies corresponding to the anticrossing of the bands, 
where a significant flattening of the Fermi surface (see insets) 
takes place. The quasiparticle energy spectra are shown in 
the insets as a function of the variable since Sj depends on 
the crystal momentum k only through C(k). 



and 712 (see dashed line in the inset of Fig. 2). Thus, 
in a wide parameter range, the hole concentration n2 
occurs near the optimum value. A typical dependence 
of the phase concentration c versus doping is shown in 
the inset of Fig. 2b. For low doping, the system is in 
a homogeneous state (c = 1). If ti > ni, the sample 
divides into droplets with two different hole concen- 
trations, ni and n2. Further increasing n, the relative 
concentration c of the phase with lower hole content ni 
decreases almost linearly, and when n > n2 this phase 
disappears and the system becomes homogeneous again. 
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FIG. 2: (Color online) (a) The density of states at the Fermi 
level versus doping in the two-band Hubbard model. Inset in 
(a): energy of the homogeneous (red solid line) and the phase 
separated (blue dashed line) states. The dependence of the 
energy in the homogeneous state, Ehom, on n has a negative 
curvature if ni < n < n2 . In this range of doping, the PS state 
becomes more favorable since its energy, Epa, is lower than 
Ehom (see inset in (a)). The values ni,2 are indicated in the 
main panel of (a) by black vertical dotted lines. The point n2 
is near the peak in the DOS. (b) The critical temperature, Tc, 
of the superconducting transition versus doping level n for the 
homogeneous (red dashed line) and PS (blue solid line) states 
calculated using Eq. The Tc{n) of the homogeneous state 
decreases fast when the doping n deviates from its optimum 
value about 712. In contrast to this, Tdn) in the PS state 
is a broad function and exhibits a plateau within the n\ < 
n < n2 range. Inset in (b) demonstrates the dependence 
of the concentration c of the phase with lower hole content 
versus n. The regions of homogeneous (HS) and PS states are 
indicated in the inset. Here we use the following parameters: 



Ed = 0.2 eV, Wp 



1 eV, Wdd ~ 0.3 eV, and Wp 



0.7 eV. 



Further analysis shows that the PS state occurs in the 
parameter range where Wpd < Wpp. If Wpd 3> Wpp, the 
two-band Hubbard Hamiltonian reduces to an effective 
single-band model, as is commonly accepted 13, l3| and 
the discussed cause for PS in the system disappears. 

The PS leads to a redistribution of the charge carriers 
and charge neutrality breaking. The structure of the PS 
state can be either ordered (e.g., checkerboard structure) 
or random (i.e., randomly distributed droplets of one 



4 



phase within a matrix of the other phase). For the ran- 
dom PS state, the size {D) of the droplet is determined 
by the competition between the Coulomb and surface en- 
ergies. Following the approach developed in Ref. [3, we 
can estimate D as (3 — 6)d at Wpd 1 eV, in 

agreement with the experimental data in [2|, where spa- 
tial variations of the DOS and superconducting gap were 
measured using STM. 

Critical temperature versus doping. — Hereafter we do 
not focus on a precise mechanism of superconductivity 
in cuprates. Instead, we intend to analyze the effect of 
the two-band structure of the Hubbard Hamiltonian ([!]) 
and PS on the critical temperature Tc of the supercon- 
ducting transition. To estimate we now use the BCS 
formula accounting for the Coulomb repulsion 

T, = Td exp [- 1/ {pjVp - i^c)] , (9) 

where To is the Debye temperature, Vp is the BCS 
electron-phonon coupling constant, — PjVc/[^ + 
PjVcl'n.{^/kBTi))] is the Coulomb pseudo-potential, and 
Vc ^ IfP ^ 5 eV are the Coulomb matrix elements. In 
our approach, the homogeneous state superconductivity 
can appear in different bands depending on doping, see 
Fig. 1: in band j = 1 for low doping and in j = — 1 
for higher doping. We stress that this is consistent with 
the superconducting behavior of YBCO where supercon- 
ductivity switches from 60 K to 90 K whith increasing 
the oxygen content. For intermediate doping, the su- 
perconducting state can simultaneously coexist in two 
bands, as in MgB2. Our model could also be used to 
interpret experiments [l9| showing the existence of PS in 
MgB2. However, a direct application of our studies for 
MgB2 should be justified since it is not clear whether the 
Hubbard-like description is applicable to this compound. 

The dependence of on doping n for the homoge- 
neous state is shown by the dashed line in Fig. 2b. For 
the optimum doping level, corresponding to the anticross- 
ing point (see insets in Fig. 1) of two hole energies, Tc 
reaches its maximum because of the peak in the DOS. 
Away from optimal doping, Tc{n) decreases fast with n. 
In the PS state, one of the phases retains the optimum 
hole concentration within a wide interval of doping lev- 
els. This results in a slower decrease (see blue solid line 
in Fig. 2b) of Tc when n deviates from the optimal value 
712. This provides a possible natural explanation of the 
observed Tc dependence on hole doping in cuprates. 

Conclusions. — We studied the two-band Hubbard 
model for superconducting cuprates and other multiband 
superconductors. Using the Green's functions technique, 
we calculated the hole band structure and the density of 
states at the Fermi level. The density of states exhibits 
a large peak at the anticrossing of hole bands. The peak 
is due to the flattening of the Fermi surface. Near this 
doping level, a transition from metal to insulator can oc- 
cur. Due to the peak in the density of states, Tc can 
significantly increase. For a certain range of the model 



parameters, a spatial phase separation of the hole car- 
riers between the two Hubbard bands was found. The 
discussed mechanism of PS is an alternative to the usual 
explanation of the PS, which is attributed to strong an- 
tiferromagnetic correlations (see, e.g., [H, Q). Our esti- 
mates of the spatial scale of PS is in a good agreement 
with important experimental results Using the BCS 
expression for Tc, also accounting for Coulomb repulsion, 
we show that Tc(n) is near its maximum within the dop- 
ing range where PS occurs. 
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